Tunneling electro-conductance of atomic Bose condensates 
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We consider interaction of an electron with a Bose condensate of atoms having electron affinity. 
Though states of the electron attached to atoms form a continuous band, tunneling through this band 
is strongly suppressed by quantum fluctuations of the condensate density. We adapt standard field 
theory methods originally developed for description of a particle propagating trough a disordered 
potential and present an exactly soluble analytical model of the process. In contrast with the 
standard description, we take into account inelastic processes associated with quantum transitions 
in the condensate. Possibilities of the experimental observation of the phenomenon are discussed. 



I. WHY ELECTRO-CONDUCTANCE OF 
ATOMIC CONDENSATES? 

Of course a simple curiosity first of all. But in order to 
provide a deeper motivation, we review the milestones in 
the history of Bose condensation from the first days till 
present time, trying to convince readers that this ques- 
tion appears as a natural step in the development of this 
fundamental field of research. 

In 1924 S. N. Bose [1] has calculated partition of the 
radiation field (Termodynamische Wahrscheinlichkcit fiir 
Strahlungsfeld) employing combinatorial formula for in- 
distinguishable elements. Two weeks later Einstein [2| 
has generalized this approach to the case of material 
particles forming an ideal one-atom gases and has no- 
ticed that it yields a decrease of the number of particles 
at a given velocity as compared to the prediction of the 
Maxwell law. Six months later, Einstein had stated di- 
rectly that Bose distribution below a certain temperature 
is incompatible with the conservation of the number of 
particle, unless one assumes that some fraction of par- 
ticles is condensed in the state with zero momentum. 
These classical results were purely thermodynamical be- 
ing based only on the notion of the size of phase volume 
per one quantum state. The Dc Broglic wave aspects of 
the condensation have not been discussed. 

The fundamental aspects of the field theory required 
for description of atomic Bose-Einstein condensates, such 
as collective quantum states [3i|, weakly interacting quasi- 
particles 1^] , and the mean-field description [g, [a, 01 1 have 
been developed in the middle of XX-th century in the 
context of superfluidity. But only at the end of the cen- 
tury it become possible to address the condensation of 
one-atom gases experimentally [3, M, Efll • 

For the radiation field, the situation turned out to be 
the opposite, - discovery of masers and lasers took place 
in years 50-th that is much earlier than understanding of 
a deep analogy between lasing from one hand side and 
phase transitions [ll| from the other hand side, includ- 
ing the Bose condensation of photons [l2| as a particular 
case of phase transitions. By now it became clear that 



the coherent laser radiation results from condensation of 
photons in the same mode of an optical resonator with 
an active media inside. But as in the case of Bose con- 
densation of gases, the first explanations of lasing were 
also almost thermodynamical being based on the rate 
equation and relations between the Einstein's kinetic co- 
efficients for spontaneous and induced radiations. The 
quantum state aspects of the radiation field have not been 
addressed. 

However, the rate equations alone cannot account for 
the coherence properties of the radiation, since they com- 
pletely ignore the important phase relations among the 
quantum states of the field. The consistent quantum- 
mechanical consideration of lasing shows that the re- 
sulting quantum state of the electrornagnetic field mode 
approaches a coherent Glauber state [13| with minimum 
ener gy-p hase uncertainty, which goes to the classical 
limit [14| when the number of quanta in the mode in- 
creases. Therefore within the kinetic approach, the state- 
ment that the induced photons are all coherent, whatever 
it means, has been taken as a sort of axiom, whereas the 
classical description 15] of the electromagnetic field above 
the lasing threshold was considered as satisfactory for all 
practical purposes. Quantum nature of the radiation 16j 
was only addressed in the context of the radiation noises 
and the photon counting statistics |17|. 

Though the theory required for the description of 
atomic Bose condensates has been developed earlier, 
from the standpoint of quantum-classical correspondence 
it is constructed in the same way as the laser the- 
ory. Based on the Ginzburg-Landau approach to phase 
transitions @, |T3|, it relies on the collective atomic am- 
plitude satisfying Gross-Pitaevskii equation. This ampli- 
tude being a quantum object by itself, still should be 
considered as a classical field in the context of quan- 
tum variables of individual atoms, - following Dirac|1.9J, 
the corresponding quantum field operators in the second 
quantization representation are replaced by a large clas- 
sical field plus relatively small quantum fluctuations. By 
making use of the Bogolyubov transformation, these fluc- 
tuations are partially taken into account and result in a 



Vlasov-like dispersion relation between the quasiparticles 
energies and the momenta, while the rest of the contri- 
bution accounting for quantum noises is usually ignored. 

The collective amplitude approach is absolutely ade- 
quate for the present experimental situation [20| , and tra- 
ditionally it is employed for predictions and explanations 
of many fascinating results observed experimentally, such 
as condensate interferen ce I2ll [2% [23j . vortices [24l [2a |. 
sound 12611 ■ excited states[23l hydrodynamic condensate 
motion [28L |29J . including that in the presence of long- 
range interatomic interactions. Moreover, this ampli- 
tude satisfying the nonlinear Schrodinger equation, in- 
vites to extend classical results of the Condensed Mat- 
ter physics on the case of atomic condensates. The 
solid-state phenomena, such as quantum wires [30|, An- 
derson localization[3l|, [S^, [33l [M, [Sa. [Sq . Josephson 
contacts |37|, nanotransistor slSSl [39l. |4Q| . metal-dielectric 
phase transition [4l|, |42, H, |4J, [4^, band structures [46l 
l47l |. etc. have already been discussed in the context of 
condensed atoms. The non-linearity of the Schrodinger 
equation enreaches the variety of possible effects that 
can be observed. Still, most of these phenomena remain 
classical from the viewpoint of the second quantization, 
since they do not address the essentially quantum field- 
theoretical aspects of the atomic Bose condensates, such 
as quantum fluctuations. 

Thermodynamical fluctuations have been considered 
at the very early stage of the development of Bose con- 
densates theorvj48l |49|, including their Dc Brogle wave 
aspects [501 . Also the kinetic aspects of Bose condensa- 
tions have already been discussed [5l|, |52, Hj, |5J|. More- 
over, following natural logics of the development of Bose 
condensate theory, one inevitably arrives at a point where 
consideration of the essentially quantum fluctuations be- 
comes crucial. One of the examples is the condensation 
dynamics. According to the Liouville theorem, the con- 
densation associated with a shrinking of the occupied 
phase volume cannot occur in a closed system, - one 
needs to invoke a dissipation and hence fluctuations as- 
sociated with the dissipation in virtue of the fluctuation- 
dissipation theorem. In the context of lasing this prob- 
lem has been intensively studied in years 60-th [55j and by 
now is well presented in textbooks [5||. Recently a con- 
sistent description of the Bose condensation dynamics in 
the spirit of the laser theory has been developed [STI. l58l| . 

In this work we consider another example, where quan- 
tum field fluctuations play the dominating role. It con- 
cerns an electron tunneling through the Bose condensate 
at the energy close to the electron affinity Ea of the con- 
densed atoms. We show that the electro-conductance 
associated with this process is indeed strongly influ- 
enced by an essentially quantum quantity - the size of 
atomic density fluctuations. The quantum localization 
mechanism [59| underlying this strong influence is well 
understood in the theory of disordered media. In our 
consideration we rely on two tools that have emerged 
from this theory: the concept of tunneling transparency 
|6(J | and the a- model technique [6l|, although both these 



tools being relevant to the case of an electron moving in 
a classical disordered potential have to be modified in or- 
der to include the case of random potentials created by 
the essentially quantum fluctuations. Indeed, according 
to the Born-Oppenheimer principle, the electron sees not 
the quantum averages of the atomic positions given by 
the mean field but their instantaneous values. The lat- 
ters, when being in generic positions, form a disordered 
media for the tunneling electron. The averaging over the 
atomic quantum state should therefore be performed only 
for the final result of consideration of the electronic part 
of problem - for the probability of tunneling through the 
disordered media |63|. 

Drawing parallels with the laser physics, one can say in 
a very broad sense, that the discreetness of atoms plays 
a role similar to the shot noise of a laser beam photons 
registered by a radiation detector. However it turns out 
that the naive analog of antibunching in photon count- 
ing, namely the binary correlation of the atomic positions 
in condensate [63], is not the parameter that governs the 
tunneling conductance. The latter depends on the size 
of collective quantum density fluctuations. Focussing at 
the simplest case, which can be exactly solved analyti- 
cally, we consider only the situation where the ensemble 
of A/" condensed atoms is in a coherent state. From the 
second quantization point of view, this state is the di- 
rect analog of the Glauber state for the laser field with 
the average number Af of photons. An analog of sponta- 
neous photons, the non-condenced fraction of the gas, is 
ignored. For soilving the electronic part of the problem, 
we also employ the simplest version of cr-model, which 
corresponds to Gaussian unitary ensembles. We believe 
that in spite of such restrictive assumptions, this consid- 
eration still remains general enough to reveal the main 
qualitative effects associated with electro-conductance of 
Bose condensates. 

We note that several questions about charged conden- 
sates have already been addressed. Transformations in- 
duced by the presence of atomic ions has been consid- 
ered in Ref . '64.] . The question of a highly efficient charge 
exchange in cold gases closely relevant to the tunneling 
conductance of condensates was the subject of Ref. [65']. 
Eigenstates and dynamics of an electrons trapped by an 
atomic Bose condensate has been considered in Ref. [6a] 
in the framework of the mean field approximation and 
the scattering length model of the atom-electron inter- 
action. All these papers were focus on the case where 
the charge is completely localized inside the condensate. 
However since the polarization energy of atoms usually 
exceed considerably the condensate chemical potential, 
by putting a charge inside a cold gas one either destroies 
completely the condensed phase or, in the most opti- 
mistic scenario, triggers a long-lasting process of relax- 
ation to a new metastable condensed configuration. In 
contrast, the tunneling conductance does not imply local- 
ization of the charge inside the condensate, and therefore 
it appears as the most gentle way of addressing the es- 
sentially quantum properties of condensates by means of 



an electrical interaction. 

The paper is organized as follows. In Sect|TT]we present 
a heuristic picture of the tunneling trough a condensate, 
which not pretending to be rigorous gives an intuitive 
idea about the process. In Sect lllll we introduce the main 
formal tools required for consideration, the Hamiltonian, 
the initial conditions, and describe a representation most 
adequate to the case of the electron energy close to Ea^ 
based on the quantum fields of atoms and a negative 
ion. In Sect llVI we consider the process with the help 
of the mean-field approximation and determine the typ- 
ical time and energy domains where the tunneling can 
occur. In Sect|V] we turn to the role of the quantum 
density fluctuations and show that they crucially mod- 
ify the mean-field picture. Trying to make this section 
as simple for reading as possible, we have moved to Ap- 
pendi JXl the most important but rather technical part 
of the calculations based on the cr-model. In Sect lVIl we 
continue by considering the regime required for the obser- 
vation of the tunneling and estimate the atomic density, 
the typical tunneling time, the typical current, and typ- 
ical losses of the condensed atoms during the tunneling 
along with there energy dependence. It turns out that 
in order to clearly observe the phenomenon one needs to 
increase the densities of the atomic condensate by a few 
orders of magnitude with respect to ones experimentally 
available at the moment. We also discuss some require- 
ments that the leads materials should satisfy in order to 
ensure efficient transport of electrons to the condensate. 
In Sect lVIll we conclude by discussing the main results 
obtained and the feasibility of experimental observation 
of the phenomenon. 



II. A HEURISTIC DESCRIPTION 

The tunneling regime implies that the wave function 
of the tunneling electron forms a wave packet of a size 
much larger then the size of the condensate, such that 
no considerable charge present at the atoms during the 
elementary act of tunneling. This is illustrated in FiglTJ 
The electron wave packet comes from the right to the 
point Ta , where the right lead approaches the conden- 
sate, and mainly scatters back. Still a small part of the 
wave packet amplitude tunnels to the left point Th, where 
the second lead approaches the condensate. The tunnel- 
ing occur through a number of the negative ion states 
corresponding to the electron location near the positions 
r„ of the condensate atoms. Each of the ionic states is 
formed as a result of the attraction of the electron by 
a short-range affinity potentials, which we chose in the 
form EaS{r—rn)- Penetration of the electron through the 
gas is a coherent process resulting from the interference 
of all possible tunnefing trajectories. 

We incorporate the well-known approach to the quan- 
tum tunneling electro-conductence trough a disordered 
media resulting from the interference of many possible 
tunneling trajectories into the quantum consideration of 




FIG. 1; Tunneling conductance of Bose condensate. A wave 
packet of the tunneling electron is much larger as compared 
to the size of the condensate, such that the electron mainly 
remains within the leads and no charge present on atoms dur- 
ing the process. The Fermi energies of the leads Ef are close 
to the electron affinity Ea of the atoms, (a) Energy diagram 
of the leads and states corresponding to various locations of 
the negative ion. (b) Wave packet of the tunneling electron. 
(c) General view of the setting. 



the condensate. Briefly, the idea of the description is fol- 
lowing. As the first step, for a fixed realization of random 
potential, one writes a formal expression 

GAi{U{r)} , 77, ra,n)GRi{Uir)} , r,, n, r,) 

in terms of advanced and retarded Green's functions for 
the probability of tunneling of an electron at energy ry 
from an initial point Tq to a final point ?■(, through the 
domain of random potential U(r), which includes inter- 
ference of all possible tunneling trajectories. As the sec- 
ond step, one takes an average over all possible realization 
of the random potential given by the functional integral 



Ga{U, ry, r„ n)GR{U, 7y, n, ra)WiU)VUir) (1) 



with a functional wight W{U) that gives the probabil- 
ity of realization of the potential U{r). This allows one 
to get rid of unimportant details related to particular 
realizations and retain only the universal contribution 
of the trajectory interference. After such averaging, the 
problem can usually be traced analytically. In one- and 
two-dimensional setting, the tunneling is considerably af- 
fected by the phenomenon of strong and weak quantum 
localization, respectively, that result from the interfer- 
ence of self intersecting trajectories. 

Tunneling of an electron through the Bose condensate 
must be very similar to the tunneling through a disor- 
dered potential 



U{r) 



AT 

E 



EaS{r ■ 



(2) 



with a random distribution of all Af atomic positions r„, 
although this process should have an important differ- 
ence, - the average has to be performed not over an 
ensemble of the all possible realization of the random 
potential, but over the quantum distribution of the con- 
densate atoms. If after the tunneling, the condensate 
remains in the initial quantum state, the analogy is com- 
plete, whereas the quantum atomic density distribution 
/3({r„}) plays the role of weight function : 



GA{{rn} , V, ra,rb)GR{{rn} , f], n, ra)p{{rn})d-^r. 



The situation gets much richer when one considers the 
tunneling accompanied by quantum transitions in the 
condensate. In this case, one has to develop an approach, 
which takes into account all possible final states of the 
condensate, including the outcomes that correspond to 
a partial or complete condensate destruction. In the 
framework of the Born-Oppenheimer separation of the 
electronic and atomic motions, the tunneling probability 
associated with the condensate transition from the initial 
state with the condensate wave function ipiniiT^n}) to the 
final state ipfniirn}) reads 

i^fn{{rn})GA{{rn},V,ra,n)r,ni{rn})d^r (3) 

X / '^}ni{fn})GH{{rn} ,r],rb,ra)i'in{{rn})d-'^r, 

and the summation over all possible final states yields 
the overall tunneling probability. 

Let us consider a heuristic picture of such processes 
comparing several situations shown in Figlll Let us as- 
sume that the condensate is initially in a pure quantum 
state given by the wave function ipmiifn})- Distribution 
of the atomic density in this state results from the super- 
position of the probability amplitudes of various atomic 
configurations, including the rare ones that correspond 
to high probability of the electron tunneling. In other 
words, for the condensate in a pure quantum state, the 
highly conducting configurations may appear as a result 
of unlikely quantum fluctuations of the atomic density. 
In Figl2] a) we show one of such unlikely configurations, 
- some of the condensate atoms form a one dimensional 
equidistant chain, which serves as a conducting trail con- 
necting the points r^ and rf,. Tunneling through such 
a chain results from the constructive interference of all 
electron trajectories going along the trail (one of them 
is shown in Figl2] d)), and the overall tmmcling time 
approximately equals to the time of tunneling between 
neighboring atoms multiplied by the number of atoms in 
the chain. 

Evidently, the probability of an exactly equidistant 
chain is vanishing small, but it increases exponentially 
when we consider configurations where the atomic posi- 
tions may deviate from the equidistant ones, as shown 
in Figl2]b). However, such deviations destroy the con- 
structive interference of the tunneling trajectories, and 
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FIG. 2; Tunneling trails for various atomic configurations re- 
sulting from quantum density fluctuations in the initial quan- 
tum state of BEG. Interference of all the trajectories that 
go (forward and back) along the same trail yields the net 
trail conduction, (a) A rare highly conducting fluctuation : 
atoms are allined in an equidistant one dimensional chain, 
(b) A more probable fluctuation : the tunneling chain is not 
equidistant. The conduction is thus suppressed by Anderson 
localization, (c) The most likely situation: many trails, each 
of which is weekly conducting, contribute to the overall con- 
ductance, (d) A trajectory can go forward and back along the 
trail. 



the chain conductance gets exponentially suppressed by 
the effect of Anderson localization. When the allowed 
non-equidistance in the chain increases, the effect of ex- 
ponentially growing probability is compensated by the 
exponential decrease of its tunneling transparency, such 
that the net contributions of different configurations be- 
comes comparable. To answer the question "which sort 
of the trail configurations gives the main contribution?" , 
one has to take into account more subtile dependences, 
and in particular the dependence on the spacial dimen- 
sionality of the problem. Moreover, one cannot exclude 
a priori the dominating contribution of numerous inter- 
fering "diffusive" trails, shown in Figl^lc). 

Dominating trajectories of different type result in dif- 
ferent degree of the condensate destruction. In fact, a 
rapid tunneling implies existence of an almost equidis- 
tant highly conducting chain. (Figl2] a). The "which 
way" question is thus answered, at least partially, and 
the condensate atoms are somehow localized on the trail 
connecting the points r^ and rh- According the standard 
prescriptions of Quantum Mechanics, the initial state of 
the condensate therefore gets reduced. The fact of a rapid 
passage of the electron through the tunneling interval can 
be considered as a sort of measurement applied to the 
condensate. Note that this measurement is indirect, and 
therefore it is not associated with a Hermitian operator of 
a physical quantity. The relevant quantity is the matrix 
element of the electron evolution operator 



fP*fn{{f7i})GB.{{'^n},t,rb,ra)A7i{{rn})d-^r, 



given by the inverse Fourier transform of the Green's 
function J dr] exp[i?7i] G/jdr^} ,r],r},, To)- As the result of 
measurement, the atoms locahzed at the tunnehng trail 
got extracted from the condensate, and the elementary 
act of tunneling thus implies evaporation of all the chain 
atoms. In contrast, in the opposite limit when the domi- 
nating contribution comes from interference of numerous 
long trajectories (Fig[2]c), the "which way" question re- 
mains without answer, and the correspondent condensate 
state reduction is much less pronounced. It can result in 
evaporation of just a few atoms per the tunneling event. 
In order to understand which scenario takes place in re- 
ality, one has to solve the problem analytically. 



III. HAMILTONIAN AND THE INITIAL 

STATES 

Let us start the analytic consideration of the conden- 
sate conductance by reviewing the relevant fundamentals 
of the quantum many-body theory of Bose liquids. The 
description rely on the quantum field theory, where the 
atomic variables are given in terms of the bosonic sec- 
ond quantization operators acting on the vacuum state 
|0) : the creation field operator $+(r) and the annihi- 
lation field operator $(r). The atomic density at the 
point r corresponds to the mean value of the operator 
n{r) — <I>+(r) $(r). With the help of Fourier transfor- 
mation one introduces the corresponding field operators 
$+(p), $(p) and n{p) = $^(p) ^(p) in the momentum 
representation. The Hamiltonian of the ensemble of cold 
atoms 



ffo 



$+(p)epl>(p)d''p + 7 / $+(r)$(r)$+(r)$(r)d''r 



involves the single particle kinetic energy given by the 
dispersion law ep and the local binary interaction char- 

^^^ given in 



acterized by a coupling constant 7 



terms of the scattering length a^ Here M is the atomic 
mass, and h is the Planck constant, which we set to unity 
hereafter unless the contrary is said explicitly. 

In the limit of large number of particles A/" — > 00 one 
can employ the semiclassical approximation for the field 
variables assuming the condensate in a coherent field 
state [13 

|A(r)) = exp |i / A{r)^+{r)+A*{r)^{r) d''r||0), 

which results from the application of the field displace- 
ment operator 



Ud {A{r)) = e'-''['^(''^*^('^)+^*('^)*('^)]'''''''. 



(4) 



to the vacuum state |0). The variational principle 
j^j^ {A(r)\ Ho - E\A{r)) = then yields the Gross- 
Pitaevskii equation 

EA{r) = €{p}A{r) + 7 \A{r)\^ A{r) 



for the displacement field eigen functions A{r). Here 
e{p) is the dispersion law tp where the momentum is re- 
placed by the momentum operator. In the presence of an 
external field the potential energy term has also to be 
included. 

The lowest energy solution of the Gross-Pitaevskii 
equation gives the fundamental state of the condensate. 
However, in zero order Born-Oppenheimer approxima- 
tion the excitations do not directly contribute to the 
condensate conductance as long as the main parts of 
atoms remain in the fundamental state. They may man- 
ifest themselves indirectly, via the mean density (n(r)) 
and the mean fluctuation of the particle number Sn = 

I (n^(r)) — {n{r)) ) . Therefore here, we ignore the 

excitations and consider the technically simplest case of 
a coherent initial state of the condensate 



exp i i / [A$+(r) + A*$(r)j d'^r i |0) , (5) 



in 



with a uniform mean density n — AA* and Gaussian lo- 
cal fluctuations Sn oc n^ " . Moreover, remaining within 
the Born-Oppenheimer approximation, we will ignore dy- 
namics of the Bose condensate during the electron tun- 
neling process, and hence the initial state of the conden- 
sate Eq. (O is the only ingredient required for the further 
consideration. 

Note that the initial state Eq.© results from the dis- 
placement of the condensate vacuum state |0) in the 
phase space of the second quantization, which is given 
by the displacement operator Ud = Ud {A{r) — const) of 
Eq.(|3]) with a spacially uniform displacement field A(r), 
whereas the vacuum state corresponds to zero average 
number of atoms and only allows for the quantum fluc- 
tuations of the atomic density. 

Let us now turn to the description of the tun- 
neling electron. The fermionic field creation ip'^{r) 
and annihilation ?A(r) operators and their counterparts 
ip'^ij)) and ipip) ill the momentum representation to- 
gether with the potential energy Eq.([2|) yield the elec- 
tronic part of the Hamiltonian 



He = / ^+b)|^V^(p)d'^p 



(6) 



Er, 



$+(r')$(r')^ (r - r') ifj+ {r)ifj{r)d'^r'd'^r 



which also includes interaction with the atoms. Here m 
is the electron mass. However, this form of the Hamilto- 
nian is not very convenient for the description of the slow 
electron tunneling at the energies close to Ea- In fact, the 
electron states relevant to the tunneling conductance cor- 
respond to the wave functions that are strongly localized 
near the positions of the condensate atoms, and therefore 
the mean values of both kinetic and the potential parts 
of Eq.® are large, and just their difference is of the or- 
der of the inverse tunneling time. In high enrgy physics 
this situation corresponds to the well-known positronium 



problem. Acting by the analogy, we introduce the second 
quantization operators of the negative ions 

vE-(r) = $(r) I ifiir - r')i?(r')d'*r' 

^+ (r) = $+ (r) / ip{r ~r')$+ {r)d'^r' (7) 



with the help of the wave function ip{r—r') of the electron 
attached to an isolated atom. In this representation the 
Hamiltonian Eq.® shifted by the electron affinity Ea 
takes the form 



He 



d'^r 



d'^riii+ir)^{r)Vir - ri)'l'(ri)$+(ri) 



(8) 

where V is the system volume. The interaction opera- 
tor V{r) has the explicit physical meaning of the elec- 
tron tunneling probability amplitude between two con- 
densate atoms at a distance r. This amplitude is given 
by the electron affinity multiplied by the overlap of elec- 
tron wave functions if centered at different atoms. 

It is convenient to work the momentum representation, 
where the field operators are related to that in the coor- 
dinate representation by the Fourier transformations 






In this representation the interaction operator 
V{p)= f d'We-'P''V{r) 



it has an explicit form 



V{p) 



K^K 



(9) 



identical for all one 
Here 



two, and three dimensional cases. 



K — ^/2mEa is a parameter in the asymptotic 
expression exp{— Kr}of the atomic electron wave func- 
tion corresponding to the electron affinity Ea, and A ^ 
2Ea/K'^ is a parameter related to the shape of the bound 
electron wave function ip(r — r„) near the neutral atom. 
Note that in the coordinate representation the kernel of 
the interaction operator reads 



V{r) = 



{2nr 



Akc"'"!''! for ID 

'■V{p)d''p={ ^Ko{K\r\)ioY2B 

AK^g-K|r| 
7r|r| 



for 3D 

(10) 
where Kq is the Bessel function. 

Numerical value of the parameters A and k entering 
the interaction Eq.Q can be found either by direct cal- 
culation in the framework of the model of (5-potential, or 
by comparison of the calculated charge exchange cross 



section with the experimentally observed values. The 
last option better takes into account the realistic poten- 
tial that binds the electron to the neutral atom. For 
Cs atoms with the electron affinity Ea = 0.39ey in the 
dimensional units one finds 



«=^5^^2.26 10nm-il 



h 



[m-] 



(11) 



The coupling constant A can be found from the value 
of the velocity dependent resonant charge transfer cross 
section 



a{v) = / 27rrsin 





Ak2 



nh^/r^ + vH^ 



-'^Vr^+^'-'t'-'dt )■ dr, 



which in the limit of high velocity v yields 

(7{V) : 



^2^2 



(12) 



(13) 



By commring this expression with the experimental data 
of Ref.feSl a = 50077^2 at the kinetic energy IKeV {v = 
3.8 lO^^i) one finds 

S > 



A == 3.9 10" 






and hence for the 3D case 



Ak2 



TT r 



-,-K|r| 



0.9 \eV\ 



exp (-2.26 lO^r [m]) 
2.26 lO^r [m] ' 



Note that the (5-potential model yields smaller tunneling 
amplitudes, which we will employ later for the pessimistic 
estimations. 

Hamiltonian Eq.® with the interaction Eq.(fTO|) and 
the initial state of the condensate Eq.® have to be 
complemented by the initial ^^(va) |0)j and the final 
(0|j^(rf,) states of the ion that correspond to the ele- 
mentary act of the conduction, where |0)j is the ionic 
vacuum state. Here we consider only the case of a sin- 
gle electron tunneling, leaving the problem of the Cooper 
pairs tunneling for future studies. Therefore the the ion 
density operator - ^(r)^+(r) and the atom density op- 
erator $(7')<l>+(r) given in terms of the field operators 
satisfy the normalization conditions 



; / ^(r)#+(r)d'^r|0)^ = 1 
Jv 

%)#+(p)-^|0h^l 

(27r) 



(A(r)l l>(r)$+(r) \A{r)) d'^r = M 



(A(p)|$(p)$+(p)|A(p)) 



d'^p 

(2^ 



= M 



(14) 
(15) 
(16) 
(17) 



IV. MEAN FIELD APPROXIMATION FOR THE 
TUNNELING DYNAMICS 



We now consider the probability of the electron tun- 
neling through the Bose condensate in the framework of 
the simplest mean field model. The results of these cal- 
culations are strictly speaking incorrect, since they get 
strongly modified when, in the next Section, we take into 
account quantum fluctuations of the condensate. Still 
the model is useful for identification of typical tunneling 
times and typical structures of the tunneling probabil- 
ity dependence on the electron energy. Though in the 
chosen regime of large electron wave packets (FigH]) the 
electron is never completely localized in the condensate, 
we will employ the dynamic formulation of the problem 
where the electron is initially localized at an atom near 
the point r^ and, as a result of sequential tunneling from 
one atom to another, reaches an atom near the point rb 
. This representation technically is much more conve- 
nient for certain calculations, while the results for the 
regime of large wave packets can be obtained by proper 
time convolutions. In other words, when considering the 
tunneling electron we neglect its Coulombic effect on the 
translational degrees of freedom of the condensate atoms. 

The initial state of the system 

\in) = <i.(r,)*+(rJe'^[^*^W+^**M]^''- 10) |0), (18) 

reveals the fact that the electron by entering the conden- 
sate in the state Eq.([5]) annihilates an atom in the point 
Ta and creates an ion at it's place. By the analogy, the 
final state of the system just before the electron leaves 
the condensate reads 

(/n| = (0|, (0|e-'/[^*^W+^**W]d'''-<l+(rb)*(r,). 

(19) 
The time dependent tunneling probability amplitude 
given by the retarding Green's function GR{t,ra,ri,) = 

{fn\ exp < —\tH > \in) has the explicit form 

Gfl,(t,r.,r,) = (0|,(0|e--^(^*"W+^**W)^'''^|.+ (r,)*(r,) 
exp I -ii / d'^rd'^rii>+{r)^{r)V{r - ri)*(ri)<l+(ri) 

|0)|0),. 



$(r<j)*+(ra)e 



+ (t \j!\A^^(r)+A''^r)Yr 



(20) 



Now we employ the representation suggested by the 
unitary displacement operator Ud of Eq.Q, which with 
the allowance for the bosonic commutation relations 



$(r),$+(r 



5{r- 



results in the transformations 



$+(/)-> $+(r')-iA* 
$(r') ^ $(r')+iA 



(21) 



of condensate field operators. This yields 
GR{t) = {Q\^{Q\[^+{n)-\A*)i>{n) (22) 

-it/drdri*+(r)(<l(r)+iA)V'(r-ri)*(ri)(<I>+(ri)-iA*) 

^$(r,)+iA)*+(r,)|0)|0),, 

where the Green's function arguments Va and r^ are omit- 
ted. 

The mean field approximation implies that one simply 
ignores the fluctuations and hence neglects the quantum 
the field operators $ and $+ acting on the vacuum state 
as compared to the average macroscopic amplitudes A* = 
A = y/n given by the condensate average density n. The 
electron tunneling amplitude Eq. (j22p in the momentum 
representation thus takes the form 



Giiit) = n 



d'^p 



rf(0|,*(p)exp[ip(r,-r6)] 



exp 



(27r) 
Ant ji>+{p)V{p)-^{p)-^ 

J (271") 



(23) 



*+(p)|0) 



Substituting the exphcit expression Eq.® for V{p) to 
Eq. (P5)) , after straightforward calculations of the matrix 
element one obtains 



GR{t,R) 



{2nY 



exp 



-it 



p^ 



\pR 



d-^p, (24) 



where R — ra — r^ and we have employed the natural 
units of the problem: the length unit 1/k ~ 0.45 lO^^m 
and the time unit (27r) h/nK ^ 5 10"'*s. The numerical 
values are given for the case of 31? setting and for the 
condensate density n ~ lO^^cm"'^. 

In the natural units, the density n becomes a small di- 
mensionless parameter ~ 10~^ that equals to the expec- 
tation number of atoms in the volume k~^ . The Green's 
function is thus proportional to n. This factor comes 
from the initial states and reflects a small probability 
amplitude to have the condensate atoms in the vicini- 
ties of the point electrodes. Evidently, this small pre- 
exponential factor increases considerably when one in- 
tegrates over the leads surfaces: for the surface of the 
order of the condensate cross-section S ~ 10~"'^^m^ the 
enhancement factor can reach ~ 10^. One may interpret 
this factor by saying that the number of tunneling chan- 
nels increases proportionally to the lead surface areas. 
Note that the transition probability scales as the con- 
densate density square, which is a natural consequence 
of the condensate coherence resulting in the constructive 
interference of different channels. 

The large parameter R allows one to evaluate the 
integral Ea.p4|) with the help of the stationary phase 
method. The saddle point is given by the conditions 



R = 



-2tpx 



(p^ + ir 



Py = Pz =0 



(25) 




FIG. 3: The tunneling probability as a function of scaled time 
t and scaled distance between the electrods R. Results of the 
saddle point aproximation on the left. Results of the exact 
calculation on the right. 
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FIG. 4: Scaled transition probability \R'^~^GR{e,R)/n\ for 
the mean field model as a function of the energy e. ID-solid 
line, 2D -dash line, 3D -dot line. 



Alternatively, this integral can be evaluated numerically, 
with the allowance for the fact that the contribution from 
the asymptotic p —^ oo corresponding to the Fourier 
transform of 6 function has to be subtracted for numeri- 
cal convergency. The dependence of |G_R,(i, R) \ on scaled 
time and distance is shown in FigOfor ID The left plot 
is done with the help of the saddle point approximation, 
whereas the right one is an exact calculations based on 
the representation of the integral in Eq. ((24)) 



GRit,R) 



2ttR 



E 

fc=0 



{-itR/2f 
k\{k-l)\ 



Ki 



{R) (26) 



as a converging series of the Bessel functions K^, {R), 



-it Ip' 



1 



in Tay- 



which results from casting of exp 

lor series followed by exact integration over p of each 
term of the series. For other dimensions just a power 
of time enters the denominator of the pre-exponential 
factor. The saddle point conditions Eg. ipS)) help to un- 
derstand why the tunneling starts at a time, which is 
Unearly proportional to R: at the real axis Px the com- 
bination —2px [Px + l) has the maximum 3v3/8, and 
hence the saddle point locates at the real momentum axis 
only when the ratio R/t is less than this number. Note 
that the contribution of only small p yields the free par- 
ticle propagator 



GR{t,R) 



(-4i7rt) 



d/2 



exp 



.R^ 

^ 4t 



(27) 



As it has already been mentioned in the beginning of 
this Section, the time dependent approach is not directly 
relevant to case of large size of the tunneling wave pack- 
ets, whereas the energy representation only has a physical 
meaning. By applying the Fourier transformation in time 
/„ dt . . . exp[— iei] to Ea. ([M)l and neglecting the singular 
terms at i? = we find the energy dependent transition 



probability amplitudes 



Gfl(£, i?)|iD = TTA/TTTexp 



2£2 V e + 1 



-R 



e+l 



(28) 



GR{e,R)\ 



2D 



GRie,R)\3D 




inRe 



which yield probabilities depicted in FigU) Note that 
in the dimensional units, the point e = corresponds 
to the electron afhnity energy, and the point e = — 1 is 
shifted with respect to this value by —nA/ {2tt) . The 
physically meaningful domain corresponds to e close to 
— 1 where the typical momentum is small as compared 
to K, which is consistent with the free particle expression 
Eq. (P7|) . Moreover, one finds a smaller cutoff momentum 
Pc ^ n^'^ for the electron by imposing the requirement 
that the phase space integral tj^ S^" d'^P equals the to- 
tal number of the electron quantum states, that is the 
number of atoms M . Going beyond this momentum is 
inconsistent with the formfactor approximation Eq.©. 
Note that the correlation effects in the interatomic dis- 



tance discussed in Ref . |63 1 



THE ROLE OF QUANTUM DENSITY 
FLUCTUATIONS 



The mean field picture gets drastically changed when 
the quantum fluctuations of the condensate density are 
taken into account. The fluctuations cause scattering of 
the electron in the course of tunneling, while the energy 
and the momentum transfer associated with this scatter- 
ing induce quantum transitions in the condensate. An- 
alytical consideration of this regime is technically rather 
involved and rely on the functional integration over the 
regular and Grassmann fields. This method, known as a- 
model, is widely employed for description of disordered 



systems, and we just adapt it to for case of averaging 
over the condensate quantum state instead of the ensem- 
ble averaging, for which it has been originally developed. 
We start with the technically more simple time- 
dependent version of Eq. ^ for the tunneling probability 



V'/n({»^n})GA({r„},T,ra,rb)'0*„({r„})d-^r 

and construct a generating function 

P(i, T,z)=J2 (»«l GAir) \fn) e'-'^ {fn\ Gfl(i) \in) , 
/" 

(29) 
where 6Af denotes the difference in the occupation num- 
bers between the initial |m)and the final \fn) states of 
the condensate. This function allows one to find several 
important physical quantities. In particular, the Fourier 
transform of the generating function 

oo 

P(,,,)./p(,,r,,)exp[-i„«-.)|<itd. 

— C30 

at z = equals to the net electron tunneling probability 
at the energy rj with the allowance of all possible changes 
of the condensate quantum state induced by the tunnel- 
ing. The derivative dzP{ri, z)\z^o gives the number of 
atoms that have left the condensate as a result of the 
tunneling. The probability of tunneling with no conden- 
sate transition corresponds to P(??, — oo). 

The explicit form of the generating function Eq. ((29l) 
corresponding to the initial state of condensate Eq.© 
reads 



P{t,T,z) 



^UU;['^ira)^+ira)e-"- (30) 






■^fr-^^^+z 



$(r6)* irb)Udc''^U^'^+in)^{n) 



-UHr 



$(ra)*+(ra)C/d|00)jO) 



where Ud is the displacement operator Eq.([4]), J\f 
J dr$+(r)$(r) is the number of particles operator. 



Ha= / drdri* (ri)<i>(ri)V^(r - ri)<l>+(r)*(r) 
Hr^ /drdri*+(ri)$(ri)t/(r-ri)$+(r)*(r) (31) 



are Hamiltonians for the advanced and retarded Green's 
functions respectively, field operators with bars corre- 
spond to the advanced Green's function, and the vacuum 
state |00)j is the direct product of the vacuum states for 

# and ^>. 

In the first line of Eq. ([50]) one recognizes the combina- 
tion (in I Ga (''■), in the last line of this equation stands 



the combination Gu{t) \in), while the middle line gives 
the expression for X^fn !/"■) 6^ {fn\ . Let us explain the 
physical meaning of the last combination in more detail. 
The initial coherent state of condensate is constructed by 
applying the displacement operator Eq.(j4]) to the con- 
densate vacuum (no atoms only vacuum fluctuations). 
If during the electron tunneling no change of this state 
occur, the inverse displacement operator restores the vac- 
uum and the operator expz/ $"'"(r)$(r)dr makes no ac- 
tion. If a change occur, the inverse displacement does 
not restore the vacuum, but yields a state \x) with the 
expectation value (a;| / $+(r)l>(r)dr \x) equal to average 
number SAf of the atoms extracted from the condensate. 
Hence, application of the operator expz J <I>+(r)<I>(r)dr 
yields the required factor e^ . 

For description of the effect of atomic density fluctua- 
tions on the electron tunneling probability, one needs to 
separate the condensate and the electron motions. There 
is a standard way to do it with the help of functional in- 
tegration. One represents the evolution operator for a 
bipartite system as a product of two commuting opera- 
tors each of which acts on only one of the parts. Each 
of these operators depends however on a common func- 
tion, a sort of external field, which serves as a variable 
for the subsequent functional integration. In this way 
one can first calculate the required matrix elements for 
the initial coherent condensate state of the atomic part 
and find as the result a weight functional for the electron 
motion. Thus obtained weight functional is an analog of 
the distribution functional W{U) in Eq.(IT]) of random 
potentials employed in the theory of disordered conduc- 
tors, and therefore further calculations can be done by 
the analogy, - with the help of a technique known as the 
non-linear cr-model. 

Details of these rather long analytic calculations are 
given in the Appendix Sect llXI The final expression for 
the generating function has an explicit form 



PiC,,z)^-2G.[mr2±^]GniRj''-''' 



/OQ pi 

dZ / dzi;e-^2(i-c-)g^(z^-^^)-i2A^C942— )^ 

(32) 



in terms of the electron energy rj 

-J. 

2 



2i^ and the frequency 

C = ^^ . Here qs = Im 1/77^ — 2 is a parameter related to 
the electron states density profile, and Green's functions 
Eg. ([32)1 given explicitly in Ea. ([84|) in fact coincide with 

that given by Eq. ipB]) where the arguments -i /^iiare sim- 
ply replaced by the combination i (77 ± iq^) /v2. Re- 
sults for the net tunneling probability V{t —> 00,77) = 
P(0,77,0) and the number 5J\f{t -^00,77) = ^^(0,77,0) 
of atoms extracted from the condensate are depicted in 
Figl5]for 2D case along with the probability suggested by 
the mean field model. As it will now be explained, the 
dimensionality of the problem plays here just a marginal 
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FIG. 5: Effect of tfie condensate density fluctuations on tlie 
tunneling probability P. The probability (solid line) changes 
with the distance R and the energy 77 according the law 
exp (— 2/?y/^) and differs from the result of the mean field 
model (dash-dot line) . The quantum state density g and the 
number 5JV oc rj~^'^'^''''^ of the atoms kicked out from the 
condensate for 2D case are shown by dot line and dash line, 
respectively. The dependences are shown not to scale, and 
the energy rj — here corresponds to the band edge 



role, and does not affect the conductance property in the 
domain, which gives the main contribution to the tun- 
neling probabihty. 

Let us give a qualitative picture of the tunneling pro- 
cess emerging from the analytical consideration. The 
mean field model considers atoms as uniformly dis- 
tributed in space and relies on the atomic density 
smeared over the space. The tunneling electron mobil- 
ity then results from the dispersion law Eq.(l9]), which 
close to the lower spectrum edge can be considered as 
the kinetic energy of a free particle with a large effective 
mass, whereas the potential energy is ignored. However 
within the Born-Oppenheimer approximation, the elec- 
tron potential energy created by frozen atoms, on the 
contrary, is a strongly irregular function composed by 6- 
like spikes placed at generic atomic positions in space. 
Were these atoms be indeed fixed in space, the tunneling 
electron would experience the Anderson localization ef- 
fect resulting from the quantum interference of many dif- 
ferent multiple scattering trajectories. At the low-energy 
edge, the electron state density spectrum would corre- 
spond to short Anderson localization lengths, and the 
tunneling would be strongly suppressed. The difference 
between this case and the case under consideration is in 
the fact that the frozen atoms are not fixed, and even 
without changing their average space positions they still 
can change quantum states by changing their momenta 
as the result of the multiple electron scattering. This 
destroys the electron trajectory interference and the re- 
lated influence of the dimensionality of the problem. It 
is therefore not surprising that the result can be found 
in the framework of zero-dimensional a-model, while the 
role of the dimensionality is reduced to purely geometric 
factors. We note that in such a situation the Green's 
functions keep almost the same form of the coordinate 



dependence just experiencing a shift of the energy scale 
and acquiring imaginary parts of the energy arguments 
corresponding to the rate of quantum transitions among 
the states of condensate. 

Suppression of the Anderson localization for an elec- 
tron tunneling through the condensate can also be viewed 
by analogy to the quantum molecular transitions. In a 
molecule with fixed atomic positions, the electronic en- 
ergies take discrete values thus forming electronic terms, 
which can be seen as a consequence of the interference 
of electron trajectories. A weak external fields is unable 
to induce purely electronic transitions without involv- 
ing the nuclear motion, unless the resonance conditions 
are satisfied, - the interference prevents electrons from 
changing of the quantum states. The resonances implies 
exact tuning of the field frequency, which for given posi- 
tions of atoms can be achieved only for a specific choice 
a pair of the electronic states. However for the poly- 
atomic molecules, when one takes into account the nu- 
clear motion with a very dense vibrational spectrum, the 
resonance conditions can be matched for many electronic 
terms at once. In other words, the quantum interference 
of the purely electronic motion gets destroyed by the vi- 
brations. The transition probability is then governed by 
the Frank-Condon factors resulting from the overlaps of 
the vibrational eigen functions corresponding to differ- 
ent electronic terms. For the condensate, the role of the 
Frank-Condon factors play the matrix elements Eq.(I3]), 
whereas the role of the vibrational energy quantum num- 
bers play the second quantization occupation numbers of 
the condensate modes. 

Let us now compare and contrast the results of the 
mean field model and the OD cr-model, illustrated in 
FiglS] Near the bottom spectrum edge the density of 
states in the mean field approximation correspond to that 
of free particle, while the allowance of the quantum fluc- 
tuations results in a Wigner semicircular density distri- 
bution. In both cases the tunneling probability is pro- 
portional to the product of the coordinate dependent ad- 
vance Ga and and retarded Gr Green's functions. How- 
ever, for the mean field model, they both are harmonic 
functions of the distance r and correspond to the free 
particle motion, while the allowance for density fluctua- 
tions makes Ga and Gr exponentially decaying with r. 
In the mean field model, Ga (^, r) and Gn {e, r) depend- 
ing on different energy variables ^ and e , respectively, 
after the inverse Fourier transformation yield two wave 
packet propagators related by the complex conjugation 
condition Ga {t,r) = GJj (i,r). The result of cr-model is 
different: Ga {rj, r) and Gr (77, r) being exponentially de- 
creasing complex conjugate functions of r both depend 
on the same energy variable 77, while the frequency depen- 
dence of the type 



(givK) faci^ors out and yields a lin- 



ear time dependence as a common cofactor of the product 
GaGr. Here g{rf) is the density of states at the energy 77 
(see Eg. ([751) ). The situation thus resembles behavior of 
the tunneling evanescent waves under a potential barrier 
in the regular formulation of this quantum-mechanical 
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FIG. 6: State densities and Green's functions for the mean 
field ID model and with allowance of quantum flucuations 
in QD (T-model. Solid line represents the state density with 
disorder, while the dashed line corresponds to the mean field 
model. These profiles can be refered to as the electron con- 
ductance bands of the condensate. 
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FIG. 8: Regular and almost equidistant spectrum of the elec- 
tron in the mean field model (to the right) corresponding to 
the plane wave eigenfunction tesults in propageting wave pa- 
chets. Irregular spectrum (to the left) with GUE statistics 
and localized eigen functions results in the evanecent wave 
that linearly increases during a finite interval of time, till the 
Heisenberg time, when it saturates. 



problem, although here the transition probability linearly 
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FIG. 7: Time dependence of the tunneling probability P with 
the allowance of the fiuctuations (solid line) and for the mean 
field model (dash line). The rate of the condensate atoms loss 
5Ni (dash-dot line). Curves not to scale. 



increases with time till it reaches the saturation value, as 
shown in Fig. [11] in the Appendixs. Allowance for the 
finite life time of the eletron in the condensate associated 
with it's returns to the leads, result in a certain trans- 
formation of this dependence shown in Figl?! although 
it still differs essentially from the tunneling decay law, 
which is usually an exponential one. 

The linear time dependence is typical of random en- 
sembles of unitary matrices, and it reveals statistical 
nature of quantum fluctuations hidden in the coherent 
states of the atomic condensate. This is illustrated in 
FiglSl Note that total tunneling probability given by the 
integral of the linear time dependence corresponds to the 



limit ^ — > in the Fourier transform 



im"(g(^)C) 



Eq.iZHD 



and always gives unity. However, the time at which the 
tunneling reaches the stationary value equals the Heisen- 
berg time given by the state density 3(77), and thus the 
smaller is the state density the shorter is the saturation 
time. This results is the consequence of the random ma- 
trix theory where all spectral properties and their time 
dependent counterparts are governed exclusively by the 
local density of states. 

One has to emphasize an important qualitative differ- 
ence between the electron moving in a disordered poten- 
tial and the electron tunneling through the condensate. 
For the condensate, one of the channels contributing the 
electron transport should be avoided, if possible. This is 
the channel when the electron makes not a virtual but a 
real transition to some point in between the leads, which 
implies the electron charge localization in this point and 
hence provokes complete destruction of the condensate. 
It does not contribute to the process of tunneling. In 
contrast, this processes does give the dominating con- 
tribution to electronic diffusion in disordered potentials, 
and, evidently, is influenced by the space dimensionality. 
In other words, the tunneling regime requires the con- 
densate size being of the order of the electron mean free 
path in the corresponding diffusion setting. Due to this 
reason, utilization of O dimensional cr-model for the de- 
scription of the tunneling looks as a reasonable approach. 



VI. THE ROLE OF ATOMIC AND THE LEAD 

MATERIAL PROPERTIES, AND THE 

TUNNELING EFFICIENCY IN NUMBERS 

In order to discuss the possibility and requirements for 
an experimental observation of the electron tunnelling 
through a condensate, we now present the typical prob- 
ability and time scales of the tunneling process in the 
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dimensional units along with the related spectral charac- 
teristics. Since the objective is to describe the electron 
transition from one lead to the other, the corresponding 
probability amplitudes should be given in terms of elec- 
tronic states, and hence the ionic point-to-point Green's 
functions Eqs. (f24l l84|) have to be projected to the lead 
electron wave functions ipi{rf,) and (pi{ra) with the help 
of Eq.Q, which relates position r of the negative ion 
with the wave function if{r — r') of the attached electron. 
The main contribution to the projection comes from the 
evanescent tails of (p; integrated over two volumes dz^dS't, 
and Az\,dSh near surface of the lead a of area Sa and near 
surface of the lead b of area Sb, respectively. 

Among all possible Lpi , the maximum overlap in the di- 
rections Za-fi -L Sa;b eventually comes from the electronic 
surface states fsiifa-b), while the overlap integrals 

playing role of the initial and the final ionic wave func- 
tions are well-localized in the vicinity of the surfaces. 
If we assume a typical size of the evanescent tails of 
the order of k~^ and take the lead areas big enough 
Sa ^ Sb ^ Kn~^, which implies many condensate atoms 
on average located at a distance < k~^ from the sur- 
face, the pre-exponential factor n in Egs. dMl [M|) can be 
replaced by a large parameter Nc — SukT^ that has a 
physical meaning of the number of tunneling channels. 
Note that "on average" is understood here not in the 
statistical but in the quantum sense. Different tunneling 
channels therefore interfere, whereas the results of such 
an interference for 3Z) are given by the integrals 

Gn^ fdSadSb'-^^^^^P^e'-^^^^-^^^^^, 

J 4:TTK\ra-rb\ 

dSadSb^^^^^l^^P^ e-'(---^K;)l-^), 

■iTTKlra - rb\ 
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corresponding to the mean field and the a models, respec- 
tively. The amplitudes cf) are now defined at the surfaces 
and are normalized by the condition J (J3* (jjdS — 1, while 
the energies e and 77 are chosen close to the correspond- 
ing condensate conducting band edges. One can move 
the big distance R = \ra — rb\ in front of the integrals 
and obtain 



GR{e,R)=iN, 



Gr{t^,R)^N, 



AttkR 



(33) 



V2 



AttkR 

where the surface functions arc now in the momentum 
representation. 

One sees that the tunneling probability depends on the 
surface electron wave functions that, in turn, are gov- 
erned by the band structure of the bulk material. For 



a broad-band semiconductor with a low electric polar- 
izability, the surface electron may form a Rydberg-like 
band spectra bounded to the surface with an energy Es 
due to the presence of the weak image force. For a small 
Eg, the wave functions corresponding to these surface 
bands are rather extended in the direction normal to the 
surface. In this case, the number of the channels involved 
in tunneling yet increases. However, it may require com- 
pensation for eventual energy mismatches between the 
band energies and the atomic affinity, which in principle 
can be done by a proper periodic external perturbation. 
One can also argue that the materials with low electric 
susceptibility have less chances to destroy the condensate 
by the Van der Waals interaction. 

Detailed consideration of the optimum choice of the 
material is a separate task deserving a special considera- 
tion. Here we will not dwell on this interesting issue and 
simply take the factors cj) in Eq. ((33|) of the order of unity. 
Still it worth mentioning that different models address 
different types of the surface wave functions. The mean 
field model by considering the tunneling electron as a free 
propagating particle of a large mass imposes the phase 
matching conditions between the plane waves of the tun- 
neling and the surface electrons. Due to the difference 
in effective masses of these electrons, the conditions are 
fulfilled near the bottoms of the relevant bands, - the 
tunneling band and the surface band corresponding to 
small momenta. The bottoms of the bands thus have to 
coincide in this case. Since now the energy only enters 
the arguments of (j), the tunneling dynamics is governed 
by the shape of the surface state density profile near zero 
momenta. In contrast, the a- model for the fiuctuations 
results in irregular localized eigenfunctions of the tunnel- 
ing band that are correlated only at the atomic scale of 
sizes, and therefore the main contribution comes from the 
surface states with the momenta of the order of k. The 
time dependence is formed by the tunneling band level 
statistics. Note that in order to be consistent within the 
tunneling regime approximation, in both cases we have 
to consider the tunneling times as short compared to the 
duration of the surface electron wave packet. 

Let us now estimate the numerical values for the 
tunneling probability and the typical tunneling time at 
the example of Gs. Fot typical dimensionless energies 
77 ^ l/R, both models yield 



P 



f N, 



yinuR 



f^nV 



KAttk'^RJ 



(34) 



that for K 
IQ-^'^m^, R 



2.26 10^ [m-ij , n 
10~^r7i amounts to 

P~2.4 10~^^. 



10"cm-3, S 



Note that the tunneling probability Eq. ([M|) is not associ- 
ated with the inverse time of the tunneling, - it is just a 
geometrical factor which gives the branching ratio of the 
small probability to make the transition to another lead 
and a much bigger probability to return back to the same 
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lead. The tunneling time corresponding to this density 
is definitely superior to the time of tunneling between 
two atoms at a distance n~^l^ , which for n ^ IO^^ctti^^ 
amounts to a completely unrealistic number 



J- tat = ^e 



10^41 sec . 



(35) 



Therefore at the condensate densities n ~ \Q^^cm~^ 
experimentally accessible at the moment, the tunneling 
electro-conductance is absolutely impossible. 

However, Tiat exponentially decreases with the inter- 
atomic distance and gets in the reasonable domain only 



for the densities 



10^ 



corresponding to the in- 



teratomic distance of lOOA, where one finds the time 



T- 



7 10~^sec 



and the correspondent tunneling probability 

F~2.4 10"^. 

Therefore in the regime of dense condensate with in- 
teratomic distance of the order of the scattering length 
a ^ bOA, the tunneling approaches a realistic limit. In 
the mean-field approximation, the typical time for tun- 
neling through all the condensate volume is of the or- 
der of uRTiat — 16msec. In order to get an idea about 
the corresponding tunneling current, one has to multi- 
ply the probability by the electron charge, and divide by 
the typical tunneling time, which again yields a small 
number J = 2.3 lO^^^A that can be increased up to 
4 10^^^ A by increasing the atomic density up to the limit 
n ~ 10^^cm~^ suggested by the scattering length. Evi- 
dently, this regime remains beyond the limits of experi- 
mental feasibility and still implies an individual electron 
counting. Unfortunaly such densities are unrealistic due 
to three-body recombination and dimers formation the 
condensate gets destroyed at a much shorter time scale 
and at much smaller densities.. 

Still the situation is not completely hopeless, since the 
electron affinity enters the exponent in Eq. (j35p . There- 
fore, the tunneling probability and the corresponding 
current increase drastically for the atoms with low elec- 
tronic affinity, such as Ca, where Ea — 20meT^|69|[70|. 
At the same atomic density, the current is of the or- 
der of J ~ 6 10~^pA, the tunneling time nRTiat — 
6 10~^^sec, and the tunneling probability P ~ 10~^, 
which can be considered as the most optimistic estima- 
tion of the tunneling process efficiency. Even at the con- 
densate density n ~ IO^^ctti"'^, one finds P ~ 10~^, 
while the tunneling time remains at the observation limit 
uRTiat — 10~^sec. Moreover, estimations show that at 
this densities the three-body recombination rate K^n^ 
with Ks'OA 10~^^cm^/sec is of the order of 0.04 sec~^ , 
which is still small enough to be ignored. However, con- 
densation of rare-earth gases is a much more challenging 
task, although the possibility to control the scattering 
length in a wide region a ~ 10^ — 10* A by making use of 



Feichbach resonances has been studied theoretically [7l| 
and addressed experimentally [72| . 

Let us now turn to the spectral properties of the tun- 
neling process, where the difference between the mean 
field and the tr-model is the most pronounced. From 
the typical size of the tunneling time one sees that the 
spectral characteristics manifest themselves in the radio 
frequency domain, in other words they correspond to the 
energy scale of 10^^ — lO^^^eV^, which practically ex- 
cludes the experimental studies based on the naive elec- 
tric circuit setting shown in Fig[T] However, it becomes 
possible to employ high-Q sources of electromagnetic ra- 
diation at the frequency domain required for the compen- 
sation of the difference between the surface band energy 
and the electron affinity of atoms and emulate in this 
way small voltage differences applied to the leads, as it is 
illustrated in Fig[9] This setting allows one to trace the 




FIG. 9: Photon-assisted tunneling of an electron from the 
surface band to the states of the negative ions followed by the 
tunneling through the condensate and the photoin assistant 
transition to the surfase state of the other lead. The frequency 
detuning ui — vi between the fields acting on the right and 
on the left leads emmulate the voltage difference IS.V . The 
Fermi energy E^f corresponds to the surface states. 



spectral dependence of the tunneling probability and thus 
to distinguish between the mean field and the cr-model 
predictions. 

According to Eqs. p3)) . the tunneling probability for 
the mean field model is energy independent as long as 
we remain near the lowest edge of the tunneling band 
and discard the effect of the surface band structure. Al- 
lowance for the density fluctuations results in an expo- 
nential square root energy dependence 



Pij].R) = 



1 / A^c 



7VV2^ \\t:kR 



e 



-2y2j)-RK 



which does not depend of the frequency C, , thus it is 
not associated with the inverse of a tunneling time and 
can therefore be interpreted as an inhomogeneous spec- 
tral line. The line width 8i] = 1/ {Rk) taken in the 
dimensional units bv ~ 5rj/Tiat strongly depends on the 
electron affinity energy and on the condensate density 
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via the dependence of Tiat on these quantities. Numeri- 
cal values of Sv estimated for Ca and Cs are depicted in 
FigfTOl which shows that the width Siy is in the domain 
of a typical diode laser lincwidth for calcium already at 
n ~ lO^^cm"'^. Note that the tunneling time ~ A/'y^ 



ln/o6v[Hz] 




ln;o6v[Hz] 



17 19 

lnion[cm] 




1 ^^ ? 

m/o/2[cm] 



FIG. 10: The inhomogeneous linewidth Sv as a function of 
the conensate density. 



suggested by cr-model (see Eq. ([78l) with Qs ~ y/rj where 
r] ^ rj — \/2) is not related with the inhomogeneous line 
width by the uncertainty relations, and in the dimen- 
sional units it is of the order of T ~ A/" {Rk)~ Tiat ■ 

We now consider the frequency dependence of the num- 
ber of atoms kicked out from the condensate as the result 
of electron tunnel. Equation (??) with the allowance of 
Eg. ((84)) for 3D Green's function can be written in the 
form 



6J\f ■ 



QsN^ 



2 {f^iqsf 
2 (AA79.)' 



GAGndV 



-d/2 



which takes into account the number N^. of the tunneling 
channels. One sees that this dependence is a rather fast 
inverse power function of the This form also shows that 
the condensate losses 6M per electron tunneling event de- 
creases with increasing number of the condensate atoms. 
By taking Qs ^ 1/-R, Nc ^ M/Rk one finds 



SAf' 



1 



Nt 



{Rk 



All 



and assuming 7 ~ 1 arrives at bN 
number is of the order of 10^ for Rn • 



jj [Rnfl^. This 
10^ and TV- 10*^. 



VII. RESULTS AND DISCUSSION 

Though our consideration seems to pose more ques- 
tions than answer, we start this section by emphasizing 



the novelty of the problem : in the most general formula- 
tion, it addresses properties of a particle that experience 
an essentially quantum multiple scattering off a material 
field in a predetermined, in our case coherent, quantum 
state. The novelty has two facets, one of which relates to 
the properties of the scattered particle that in the pres- 
ence of inelastic channels manifest some deviations from 
a similar problem of quantum transport in disordered 
media, while the other facet relates with possibility to 
reveal certain properties of an atomic condensate, such 
as the size of fluctuations, by incorporating its quantum 
field to electric circuits. 

The main result of the paper is following. In spite of 
the fact that in Bose condensate all atoms are in the same 
quantum state with a wave function extending at micro- 
scopic distances, the tunnehng electron feels the discrete 
nature of individual atoms associated with the quantum 
density fluctuations, and as a consequence, the tunneling 
gets strongly inhibited by the quantum localization effect. 
In other words, a simple mean- field picture that consid- 
ers atomic condensate as a uniform support for the mov- 
ing electron cannot correspond to reality. Though the 
allowance for irreversible transitions in the condensate 
partially destroys the electron interference and thereby 
facilitate the tunneling, it does not remove the localiza- 
tion effect completely. In the accordance with this pic- 
ture is the fact that the state density profile of the tun- 
neling electron band has the Wigner semicircular form 
in contrast to one corresponding to a free particle prop- 
agation. Superposition of this profile with the profile of 
the localization length determines the spectral line of the 
conductance efficiency. 

The main message emerging from the consideration of 
relevant numbers is that the electron tunneling through 
a condensate is not a completely impossible process, and 
it can, at least in principle, be realized in a regime that 
approaches the domain of parameters currently available 
experimentally. The main problem is to identify atomic 
species that have weakly bounded negative ions. The 
typically large (~ 10'^ A^) charge-exchange cross-section 
of negative ions in normal conditions that can yet in- 
crease for cold collisions are still not sufficient to ensure 
the efficient electron tunneling transport at typical inter- 
atomic distances of currently existing condensates. Nev- 
ertheless, since the charge exchange rate exponentially 
increases with decreasing electron binding energy in neg- 
ative ions, in some cases of closed shells atoms with low 
polarizability, like Ca for instance, the tunneling regime 
can be achieved in principle, and the question is thus get 
moved to the possibility to condense these atoms. 

Note that under condition of the exponential depen- 
dence of the tunneling probability on the interatomic dis- 
tance, the role of careful calculation of specific atomic pa- 
rameters, such as the actual shape of wave function of the 
attached electron, based on standard methods of atomic 
physics may becomes important. It may change essen- 
tially the estimated required atomic density, as it was 
the case for the presented calculations performed with 
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and without the parameter A. One should also have in 
mind the case of He atoms with closed s-shell that do 
not have electronic affinity at all since their polarization 
potential is not strong enough to accommodate a bound 
state of electron. Still they are able to create a disordered 
potential for the tunneling electron at the energies close 
to zero. This potential will also modify the state density 
profile at this region. All the calculations based on a- 
model thus remain valid, and the only difference will be 
in the sizes of the time and energy units, that presum- 
ably should correspond to a more efficient tunneling. Of 
course, one needs to know precise shape of the electron 
binding potential of atoms and the corresponding scat- 
tering length in order to make more accurate predictions. 

Not a less important question is the efficiency of in- 
terface between the electric leads and the condensate. In 
order to ensure an efficient tunneling from the leads to the 
condensate without serious destruction of the latter, one 
needs to have electronic surface states rather extended 
in the direction perpendicular to the surface. Reliable 
determination of such states is a subject of serious nu- 
merical work, but at the first glance one can hope that 
the broad-band semiconductors may offer such a possi- 
bility. One also can consider the possibility to employ 
carbon nanotubes for the purpose, having in mind that 
the questions about the dielectric properties Ref . [TJ] and 
Casimir potentials for the atom-tube interaction Ref. |74| 
have already been considered in details for such systems. 
However, for the case of nanotubes one has to expect con- 
siderable decrease of the number of the tunneling chan- 
nels, which is a direct consequence of the small size of 
the surface area. 

It also worth mentioning that application of periodic 
external fields may not only compensate for existing mis- 
matches among the surface state band and the conden- 
sate tunneling electron band, but also offer additional en- 
hancement of the transition between these bands when 



goes through highly excited and therefore highly ex- 
tended intermediate surface states. Moreover for a mul- 
tifrequency resonance setting, it can also serve as an ef- 
ficient tool, which reduces the influence of the thermal 
effects in the leads by satisfying strict resonance condi- 
tions at the final or intermediate states of the excitation 
cascade. 

Looking at prospective, one can believe that the 
present consideration may also become relevant in the 
context of condensed Rydberg atoms, where the conduc- 
tance might exist due to the partially occupied bands of 
delocalized Rydberg states. However, in order to consider 
this case, one needs to allow for the electronic many-body 
phenomena, such like Cooper pairs, that have been omit- 
ted in the present work. 

One should also have in mind that there exist other 
physical systems, such as Bose condensates of excitons at 
the interface between two semiconductors, where similar 
phenomena can take place. 
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IX. APPENDIX 

Here we show how to separate the atomic field op- 
erators from the ionic ones. Prior to the separation, 
we make some transformations that simplify the calcula- 
tions. Consider Taylor expansions of the time dependent 
exponents in Eq. (|30[) . Each term of the expansion is a 
sequence of the operator products 



J 



k\ 



nd. 



, *+(r3)$(r3)V^(r3 - rz) $+(r2)*(r2)'l'+(r2)$(r2) V{r2 - ri)$+(ri)\E'(ri) 



(36) 



r 



where the left argument of factor must coincide with the 
right argument of the next operator, as it is emphasized 
in Eq. (|36p for the second argument by the underlining 
brace. This condition is the direct consequence of the as- 
sumption that we have only one tunneling electron. The 
condensate operators also correspond to the same coordi- 
nates, and hence they can be grouped together when the 



Taylor series are summed back to the exponent. In other 
words, operators <^^[r) can be moved to '^{r) comprising 
the number of particles operator in each point where the 

ion field operator ^ is present. The same is valid for ^ 
operators. Equation ([30]) thus takes the form 
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gi/(A*+(r)+A*|.(r))drgz/dr*+(r)<l(r)g-i/(Ai+(r)+A**(r))dr(|,+ / ^N,!/ ^N /gy-j 

I 



The general idea underlying the separation of the ionic 
and condensate fields rely on functional integrals, and in 
particular, on the functional Fourier transformation. The 
main steps of this procedure are as follows. A functional 

il I X{x)Y{x) ) depending on a product of two commut- 



ing operator fields X{x) and Y{x) can be written as a 
functional integral 



Vu{x)n (x{x)u{x)\ 6f {y{x) - u{x) 



where the functional Dirac (5-function has the Fourier rep- 
resentation 



Here u{x) and u*{x) are two independent real- valued 
functional variables. By combining these two expressions 

one represents fl (X{x)Y{x)] as a functional integral of 
a product, where one factor depends only on the operator 
field X{x) and the other only on the operator field Y{x) 



\Y {x) — u(x)\u* {x) 



rj(x(a;)M(x))e'2-/dx 
In particular, two factors in Eq. (|37] 



Vu{x)Vu*{x) 



ir / drdri* (ri)V(i — ri)>I'(r)<l+ (r)<i>(r) 

and 

— it/ drdri'I'+(ri)y(r — ri)*(r)'l+(r)<i>(r) 

can be represented in the form 

ir/drdri* {ri)V{r—ri)4'(r)u(r) 
gi2^/dr[4+(r)4(,.)-u(,.)]u-(r)p^(^P^*(^) 

and 

/—it/ drdri>I'+(ri)V(? — ri)%(r)v(r) 
gi2^ / dr[4+ {rmr)-v{r)]v' {r)-py^^)-py* (3.)^ 

respectively. 



By regrouping the operators, one obtains Eq. (|57|) as a 
functional integral of a product of three matrix elements 
each of which is taken by the corresponding vacuum state 



P(i,r, 2)= / Vv{x)Vv* {x)Vu{x)Vu* {x) 



where 
Gr — 



(38) 



g-i2^«(r)«*(r)g-i2^«(r)«-(r) ^^^^ q^ (^) q^ (^) 



\]/('^^~)gi-r/d'-d'-i'I' (ri)V(r-ri)i'(r)u(r)^ (rf,) 10) 

<I'(rfc)e^'*''''*'"'*''^*^^'"^'^*''^''^^*'''^"'''^'I'+(ra) 10) 



are the advanced and retarded Green's functions calcu- 
lated over the vacuum states of the ion, and the matrix 
element over the condensate vacuum 

{Cd) = (0|e-'-^('^*^W+^**W)'i'^$+(ra)<l(r,) 

i2iT-/dr#+(r)*(r)M*(r) i/(A'l+(r)-|-A**(r))dr 
g2/dr|.+ (r)4(r)g-i/(Al>+(r)+A'*(r))dr|,+ / ^-j^/ N 
i2Tr/dr<l+(r)<l(r))j*(r) i/(A<l+(r)+A*<l(r))dr \q\ 

The bar over the Green's functions Ga and Gr shows 
that they are calculated for the operators that include 
fields u(r) and v(r), respectively. 

We calculate (Cd) with the help of two types of oper- 
ator relations : the Baker-Campbell-Hausdorff rules 



^A+B ^ ^B^A^[A,B]/2 

e\^ = e^e^eI-4'^1 



(39) 



for the operators with c-number commutator, and the 
relations 

g2/dr|.+ (r)*(r)^+^^>||,^^>|g-z/dr*+(r)*(r) ^ e^l>+(r) 
g2/dr<l+(r)*(r)|)+^^-j|,j-^-j^-z/dr*+(r)|.(r) ^ e"^l>(r) 

(40) 

for the bosonic operators. With the help of Eq. ((40|l one 
can move the factors depending on the number operators 
towards the right vacuum state, where they just result in 
factors 1, and obtain 

(0|e-'/(^^"'""'"^'*+«+^*^'""*"''*W)^'-4+(r,)$(r,) 

i/(A$+(r)+A*<l(r))dr -i/(Ac^l>+(r)+A*e"^<l(r))dr 



^+{n)^{n)e'^^^'' 



+ i2^^*(^)S,+ 



*+(r)+A*c 



-^-i2,ri,*(r; 



<i>(r))dr 



|0) 
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Ignoring the quantum fluctuations in the pre-exponential 
factors yields 

Aj(Ai+(r)+A*i(r))dr ■~iJ(Ac'i+(r)+A*c-'^(r))dr 
i/(Ae'+'2''''*''''4+(r)+yl*c-^-'^""'<'''<i>(r))dr |q\ 



With the help of Eg. ([39)) the exponential terms can be 
normally ordered, that is rearranged in such a way that 
all $ are to the right with respect to $"''. All the terms 
containing normally ordered field operators vanish. The 
c-number commutators appearing in the course of the 
normal ordering yield a factor e™^" ^''■*''" '■'"-'-', which reads 

w{u* {r),v*{r)) = n f dr (e'^^"* (''^ + e'^'^^* ^''^ - 2+ 

e^(l-e'2'^"*M)(l-e'2'^"*M)). 
The functional Fourier transform 

Wu,v = /"pv*(a;)I?w*(x)e"'-2'^'-'''^''["('')"'('')+('')'''('')l, 

(41) 
of e*"^" ('')'" ('■') plays role of the weight factor for func- 
tional integration over the fields u{r) and v{r). However, 
it has a structure inconvenient for calculations. Techni- 
cally is much simpler to perform the functional integra- 
tion over the fields u and v first, and than integrate over 
u*{r) and v*{r). 

Equation ([38| thus takes the form 

P(t,.,=) = „^/l,.-(.)P.»I,..(.)P.(x, 

Ga (t) Gr u\ q^-'^^'^ f '^^M^)'^' i^')+i^')^' i^)] 

One can perform the Fourier transformation over times t 
and T and write 

P(e, £,, z) = n^ f Wu,vT>v{x)Vu{x) (42) 

where the Green's functions can be given explicitly 

Ga (0 = (0| ^(ra) (e - HuY'^^in) |0> (43) 

Gr (e) = (0| ^(rfc) (e - i?,,)"' ^+(ra) |0) 
in terms of the field operators 

Hu = / drdri* {ri)V{r - ri)i'{r)u{r) 
Hy = / drdri*+(ri)F(r - ri)^(r)u(r). 



that differ from the Hamiltonians Eq. (|3T|) having classical 
fields u(r) and v{r) instead of the quantum condensate 
fields $+(r) and l>(r) . _ _ 

In order to find the average of Ga (0 Gr (e), we follow 
the prescription well elaborated in the theory of disor- 
dered systems and express the Green's functions with the 
help of functional superintegrals over supervector fields 
A(r) and nir). For the advanced and retarded Green's 
functions these fields are grouped to the line vectors 

A(r) = (A,(r) A/ (r) ) ; ;. (r) = ( /ifc (r) fif{r)) 

respectively, and the conjugated columns vectors 

This integral over regular (bosonic) functional vari- 
ables Xb(r) and /ifc(r) and anticommuting grassmannian 
(fermionic) variables Ay(r) and fJ-f{r) read 



GR{e) = J VX{r)VX{r)e'''^Xf{ra)Xfin) 

Ga (0 = / I?A*(r)I?M(r)e-'^« V/(^6)M/(^a), (44) 



where the conjugated variables denoted by bars and the 
actions have the form 



/ drA(r)eA(r) — / drdriX(ri)V{r ~ ri)X{r)v{r) 



Ss 



5*5 = / dr/i(r)^/i(r) — / drdriii{ri)V{r — ri)ii{r)u{r). 

(45) 

Here the products A(r)A(r) and ^{r)jl{r) stand for the 
scalar products Xb{r)Xb{r) + Xf{r)Xf{r) and ^b(r)/I{,(r)-|- 
fif{r)'JIj:{r), respectively. The functional differentials 
VX{r)VX{r) stand for VXbir)VXb{r)VXf {r)VXf{r), and 
'D^{r)'D'Jl{r) for 'Dfib{r)'D'JIij{r)'D^f{r)'D'Jlj:{r), and the 
operator V reads 



^y^^ "^^l-i Vir-n) 



We substitute Eq.lgl]) with the actions Eq.fgS]) to 
Eq.((321) and perform the integration over the fields u 
and V first. This integration yields Sp functions with 
the arguments 27rw*(a;) — J driX{ri)V{r — ri)A(r) and 
2t:u*{x) + J drifi(ri)V{r — ri)]l{r). Performing the inte- 
gration over the variables u* and v* one finds 
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(46) 



P{^,e,z) =n^ / m(r)PA(r)e'/'i'-^W^^«A/(r,)A(rb) / Vfi{r)VTi{r)c-'f^-^^^^'>^''^^'>^if{n)Jlf{ra) 

I 



Thus far the calculations were exact. At this point density n, which is a small dinicnsionless parameter. By 

one has to make an approximation, since exact calcu- the second order Taylor expansion of the terms entering 

lation of the functional integrals are further impossible. the exponents in the second and third lines of Eq.(|46p we 

We notice that one fold integrals with the exponentially obtain 



decreasing kernel V{r~ri) are of the order on the atomic 



Pi^,e,z) =n^ / m(r)m(r)I?A^(r)I?/i(r)A/(ra)A(rfc)Ai/(rb)7I/(rJe'/'^''^«^^W-'^'i'-^W«^M 

in J'[A(ri)y(r — ri)A(r) — /j(ri)V(r— ri)7i(r)ldridr+c^n J dridr2dr[\{ri)V (r—ri)\(r)fi{r2)V{r—r2)jL(r)] 
g— -f /dridr2dr[A(ri)V'(r— ri)A(r)A(r2)V'() — r2) A(r)+/j(ri)y (i — r i)JI(r) fj.{r2}V (r — r2)JT{r)] 

I 



Let us rewrite this expression in a more compact way by 
introducing two supervectors in the functional coordinate 
space - the column vector | (^) defined by their coordinate 
representations as 



(47) 



<^(r) = (r|</.) = 



A/(r) 



and the corresponding line vector 



<^(r') = (0|r') = (Ab(r') /i, (r') A/ (r') ^,f{r')) 



In this notations Eq. (|T7)) takes the form 



P{te,z) ^n^ / V\^)V{(f>\ Xf{ra)X{rb)^lfirb)Tifira)e^p [{(f>\ (iC + i?/S,+inS,V 
exp 



exp 



--(1 - e^) / dndradr (0 (ra) V{r2 ~ r)ll,4> (r)) {cp (n) V{ri ~ r)S,(/. (r)) 



--(1 + e") / dridradr (0 (ra) V{r - ra)^/. (r)) (0 (n) V{r - n)(j) (r)) 



where two matrices 



10 



{r'\V\r)=V{r'-r)\ [J J J [J | and (/| S, |r) = (5(r' - r)E, = J(r' - r) 
1 



/ 1 
0-100 
10 

Vo 0-1 



(48) 



and two variables rj = ^i^ and ( — ^2^ have been in- intcgrationby all eight components of the vectors | 
troduced. Here by V \(j)) V {^(j)] we denote the functional \4'\- 



and 



19 



Let us rewrite Eg. pS)) in an even shorter way by in- 
troducing an object similar to the density matrix - the 
direct product 

P=\'4>)® {(j)\ . 

This allows one to put Eq. ((47| in the form 



P(C,?7,z) =n^ jv\4>)V{cp\ \f{ra)\{n)tif{n)jlf{ra) 

gSTr(iC+ir;S,+mS,V)p /^g-j 

g-f (l-c")STrpVS^,pVS,-f (l+c")STrpVpV 

The symbol STr stands for the regular trace operator 
with respect to the coordinate dependence and with re- 
spect to regular vector components and implies change 
of the sign for the Grassmann components. By Str we 
denote the supertrace operation, which involves only the 
matrix variables and does not imply integration over the 

I 



coordinates, it amounts to the difference of the traces of 
the boson-boson and fermion-fermion blocks. 

The next step is application of the Hubbard- 
Stratanovich separation of the terms quadratic in p. This 
can be done by introducing an superintegral over a su- 
permatrix Q which, generally speaking depends only the 
coordinate r, that is (r'| Q \r) = d{r' — r)Q with 



Q 



Q> 



Q = 



AfcAb l^AbPi, ^XbXf '^XbfJ.f 

V^bAfc ^fibl-^b ^fib^f ^t^bP-f 

^\fXb VAj/^b ^\f\j ^\jfj,j 

WlJ.f\b '^M/Mb ^t^f^f ^M/M/ 



(50) 



where the matrix elements Qi^ — Qij{r) are the regular 
functions of r in the diagonal 2x2 blocks (boson-boson 
and fermion-fermion) and the Grassmann functions of 
r in the off-diagonal (boson- fermion) blocks. The third 
line of Eq. ((49|) then takes the form 



exp|-JSTr[(l-e^)pVS^pVS, + (l-he^)pVpV]} = f VQ{r)exp 
where (r'| Q' \r) = 6{r' - r)Q' with 



Q' = 



I -Q\b\b 

I - -/2 



ie^/'QA.p. 



-Q 

zl2. 



AbAf 



le 



■.12 



--STrQ^-inSTrVQV 



(51) 



Q\b^ls \ 



le Qfj.b\b Qp-bfj-b IS Qp,bXf Qfj.bfj.f 
-QxfXb ie"/'QA,p, -Qa,a, ie^/'QA, . 

I 



In these notations Eq. (f49|) reads 
P{C, V,z)^n^ ( V \J) V (01 1?Q(r)e~^^^'^' 



V(r„)A(r,)^H^^)M/(^a)e^^-('';+'"^^+'"^^^-'"^«')^. 

(52) 

We note another possibility, - the factors e^/^ can be 
moved from the matrix Q' to the matrix Q . This trans- 
formation has a unit Berezian (the supermatrix analog of 
Jacobian) since factors e~^/^ and e^/^ arising from four 
bosonic and four fermionic variables cancel each other, 
and Eq. ([5^ takes the form 



P(C, v,z)^n^ J V\^)V (01 VQ{r)Xf{ra)Xf{n) (53) 
^/(r,)^^(r,)eST'-('^+"'^^+'"^^^+^"^^^^^^)'' 



exp 



-STr (1 



.-^/2 



)Q+(l-e-"/2)S,QS, 



where we have employed the fact that in this represen- 



tation Q-> i(S^Q + QS^ 



-z/2 



(Q - SzQS^) and 



The main features of the tunneling conductance can 
be found in the framework of the simplest version of 
this technique, which ignores the coordinate dependence 
Qi,j{f) — Qi,j of the matrix elements Eq . ([50)1 . .7a| It 
is applicable when the perturbation of the mean field 
regime by the condensate density fluctuations is strong 
enough, in complete analogy with the situation known in 
Quantum Mechanics, where a strong time-dependent per- 
turbation induce quantum transitions between slightly 
detuned levels as if it were resonant. In such a situ- 
ation, only the integral of the time dependence enters 
the result, while the particular form of the time depen- 
dence plays no role and can be taken uniform. Here, 
by the analogy, we have to compere the interactions re- 
sulting in two last terms in the exponent in the second 
line of Eq. ([5^ . Since we consider tunneling at a long 
distance i?, the typical difference of the momenta of the 
advanced and retarded Green's functions is small, and 
according to Eq.® the typical variation of the interac- 
tion associated with the term STrnSz V can be estimated 
asNdpV(jp)R^^^AfKRr^. The fluctuation-induced cou- 
pling associated with the term STmVQ can be esti- 
mated from the condition ^STrQ^ - ^SivQ'^ ~ 1, that 
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is Q ^ 1/y/N, and hence this coupling is of the order of 
A\/7\f. Therefore, the fluctuations can be considered as 
uniform for J\fR~^ ~ nR <C 1, which we assume to be 
the case for the system under consideration. 

The assumption of spacial uniform matrix elements al- 
lows one to employ the momentum representation where 
two momenta p and q corresponding to the advanced and 
retarded Green's functions, respectively. This simplifies 
considerably the integral over I?Q(r) in Eg. ([5^ . which 
now becomes not a functional, but just an 8-fold integral 
over the regular and an 8-fold integral over the Grass- 
mann variables 

dQe~^STr((i+c--/^)Q+(i-c--/^)£,Q£,)' (54^ 

^/b)-^/(p)A*/ ('7)7^/(9) 

The pre-exponential factors can be written as derivatives 



i nV{p) 



Ja 
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(55) 
(56) 



and 



l^f{<l)f^fi<l) 



9.h^0 j Ij^nVip') 4>{p')JaHp'). 



Jb = Jbjb 



i nV{q) 
(27r)' 



Jb 



V 







1 



S{q^p'), 



(57) 
(58) 



which allows one to write Eq. (|54p in the form 



p(C,7?,z) = al,7, 



dpd<7 ,j^(„_ 



^° ' (2^)^"^^ 



R(p-q) 



V(I}V(J3 

V{p)V{q) 

(59) 



/dQe-TSTr((i+c-^/=)Q-(-(i-o-^/=^)g.QS.)' 

/ ^^*(p)(C+i'?S.+n§.1^(p)+"V(p)(Q'-J„ + Ji,))0(p) 



One can substitute now explicit expression Eq.Q for 
the interaction potential, (we employ the scaled units 
•q — > TynA, C — > C,nK,t — > t/nK for t, 77 and C,, re- 
spectively, that have been introduced in Eg. iPH) '). and 



consider domain 77 ~ —1 by making the transformation 
77 —> 77 — 1, along with trhe assumption 77 ^ 1 implied by 
the physical requirement p, q << k. This yields 



PiCv,z)^dlj,^o 



dpdg ^,nip-,) I p^p^ 



s2d' 



dQe" 



i2n) 

-=^STr((l-|-c-^/^)Q+(l-c-^/^)£,Qg,) 



g /^</'(p')(iC+i('?-p")s.-i {Q'-Ja + J,))4>{p'} 

where ^ and 77 were neglected when compared to unity. 
After having shifted the matrix Q by C^z? one obtains 
a form of this expression 



p{C,v,z)^dlj^ 



dpdq ^i 

(2^) 



2rf 



R{p-g) 



V(j)V(j) (60) 



dge" 



;(4CStr£,Q+(l+c-")StrQ^ + (l-c-")Str(Q£,)^) 



p / H!y^'^(p')(i('?-p'')S.-i(Q'-.7a+J0)'«p') 

convenient for further calculations. 

We did not yet specified the precise form of the ma- 
trix Q that allows one to change order of the integra- 
tion in Eg. ([SO)) and evaluate the integral over VcjjVcj) first. 
The standard approach suggests to express this matrix 
in terms of the generators of the superalgebra osp{l, 1|2) 
and make use of the parametrization 



g = TPT- 

where P is a supermatrix 



P = 



/ qi Ki 

92 7C2 

Ki iqi 

\ K2 iq2 



(61) 



(62) 



and r is a pseudounitary rotation 

f = UM 
with 



/ cosh Q e'* sinh 8 



M = 



i* 








sinh 8 cosh 8 

cose* ie^'^sin^ 

ie"''^sin6' cos 6* 

f l-^ 0_ -a \ 

1 + '^ o_ -ip 

a 1 + ^ 0_ 

V i/3 1 - ^ / 

In this representation the integration measure (Haar 
measure for the chosen parametrization) reads 



U = 



dQ = dPdT 



dT = dcosh2e d cos 261 



(63) 



d^dip dad(3dadj3 



47r2 (cosh 28 -cos26l)^ 
dP = d(7idg2dgidg2dKidK2dKidK2- 
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We apply this parametrization to the generating func- 
tion in the form Eq. (|53p and perform the transformations 
Eqs. (|54ll59|) . Evaluation of the integral over VcjfDcj) yields 



P(C,?7,^) = 5}„„/,^o 



dpdg ,j^,„_ 



(2^) 



2d 



R{p-q) 



AQe 



^N 



-StrQ^+ ^~° ' str(QS^)^-CStrS^C 



- I i0^Str\n[n-p'-'-Q' + {j^^J^)^,] 



By substituting Q' = — -wE^QyEz and by moving y S 
under the supertrace operators one finds 



P(C,77,^)=5}^,j,-.o 



dpd? i^/p_ 



(2^) 



2d 



i?(p-g) 



-7^ i±£^StrQ2 + ^ 



Str 



dQe 

^ i7i^Strln[,,-p'^-Q+(j„-J,)] 



(Qg,)^-CStrS,Q) 



(64) 



where by the analogy to Ea. ((28|) the Green's functions 
(depending on the matrix Q) read 



Ga,r (q 

Gam [Q 
Ga,r (q) 



Wv-Q 

_ nKo{±iRy/,j-Q) 



2tt 






for ID 
for 2D 
for 3D, 



and the choice of signs has to be done such that to assure 
the convergency. 



which after taking the derivatives 

-.2 /• _ r Vdp^ 



p(c,.,.)^-ydQe-^^-'""-'^-«i 

A/'CStr£^Q-Af ^+°~° StrQ''-A^ ^-°~" Str(Q£^)^ 



dpdg 
(2^0^' 



,iR{p-q) 



1 



1 






allows one to explicitly perform the integration over dpdq 
and results in 

P(C,r;,z)= UQe-^^""'''''^^-^''-^^ (65) 

gAACStrS,Q-AAi±|^StrQ^-AAi^^§^Str(QS,)^ 

Strj,G(Q) StrJfcG(Q) , 



By moving T under the supertraces and superdetermi- 
nants, one can obtain Ea. (|65p in the form 



P(^, 77, z) = / dPdPe ■' 2(i^ ^ ' -^ J 

g-AA§Str£'^P-AA ^+°~' strP^-AA ^~°~^ Str(p£;)^ 

Strj,G(p) Strj6G(p) (66) 

convenient for integration over the matrix P . Here the 
integration over dPdT has to be done according to the 
measure Ea. (|63p . the matrices with primes read 



S' 



/ cosh 20 -e2"i'sinh2e 



-2i* 



sinh2e -cosh 29 







cos 26* -ie-2'"*sin26l 
ie^''^sin26' -cos 26* 



Ja = 



-^e2**sinh2e 
—aa sinh 



Jb 



V 

/ aa cosh^ 8 
^e-2i*sinh2e 
— acoshGcos^ a sinh 0e^'* cos 

V -aie^'"^ cosh 6 sin 9 aie2'(*+"^) sinh 6 sin 

^/Scosh^e -^e^'* sinh 26 -;5e2'(*+'^) sinhesin^ -;5ie2^* sinh 9 cos 61 \ 

-'PPsYoh^Q -^e^'-^ cosh e sin 6* -;9i cosh 9 cos 6* 



— a cosh 9 cos aie ^"*' cosh 9 sin 6* \ 

-a sinh 9e-2ii> cos 9 me-2'(*+'^) sinh 9 sin 9 
(1 + aa) cos^ 9 
ie 



-le 



.2i01+aagj^26l 



2i01ipgj^26l 

(1 + aa) siv? 9 



( 



^e-2'*sinh29 



V 



_^g-2i($+0) gjj^j^ Q gi^ g, ^e-2i0 (,Qgj^ Q gjj^ gi 



/3ie-^'* sinh 9 cos 61 



-/3icosh9cos^ 



(1 - /3/3) 

;„2i0 1-/3/3 



sm 



le 



-2101-/3/3 



sin 26* 



sin 26* 



(1 - /3/3) 



cos 



/ 



(67) 



(68) 



(69) 



22 



and the terms in the exponent have the exphcit form 

Strt'^P = {qi ~ 92) cosh 29 - i (gi - §2) cos 26* 
Strp2 = qf + q1+qf+q^ + 2kiKi + 2k2K2 

Str (pK) = {li + <&) cosh^ 26 - 2qiq2 sinh^ 29 + (g? + q^) cos^ 261 + 2qiq2 sm^ 261+ 

2 (kiKi + K2K2)cosh2ecos26l - 2i ('c^'^'^+'^^ksKi + e~'^'^'^"^'l''^KiK2\ sinh29sin26l. 

I 



(70) 



The integral over dP can be evaluated by the saddle 
point method. The main contribution, generally speak- 
ing, comes from a saddle point where the matrix ele- 
ments qi,q2, qi,q2 of P are shifted to the complex plane, 
whence the matrix P has a diagonal form Ps no longer 
corresponding to Eq. (j62p and independent on Grassmann 
variables. The latter, as well as the deviations of the 
regular variables are considered as perturbations. Note 
that the pre-expnential factor does not change as the 
result of integration over these perturbations, since the 
Grassmann integration around Ps compensate for the 
pre-factor resulting from the integration over the regu- 
lar deviations. Also note that the among different pos- 
sibilities, the saddle point has to be chosen in such a 
way that the integral over dT converges. Moreover, the 
choice qi = — 92 = i^i = ^ i92 = i^s at the saddle point 
allows one to get rid of the large terms MStiPg and 
STr ln[77 - p''^-P's], and Eq.® takes the form 

P(C, ry, ^) = n^ / df e-2-^«'«= (=°'''^ ^""-'"^ ^e) 



StrjaGAM {Ps) StrJbGA.R [Ps) 



(71) 



where the choice between the advanced and retarded 
Green's functions depends on the position of the matrix 
element qi = iig^ on the complex plane. 

Since Strj^G = aa(cosh26 - cos 261) '^""'^^ - 
sin^ 6l(Gfl - Ga) and Strj^G = /3^(cosh29 - 
cos 20) '^^"'^^ + sin^ 6'(Gfl - Ga) by integration over 
d^d(j)dad(3dad(3 , one can explicitly find the pre- 
exponential factors in Eq. ([7T|) 



Strj^G Stride 



TT^ 



-(cosh 26 - cos20f(GR - Ga? 



with the Green's functions 



/o , nc V , I 4i ^ 
^A ^ 2V^+ig. ^« 
^ ^ nKa{ — iR\/ ri+'iqs) 
(jA -^ 



.iIiVr]-qs 



2V»?-i9s 



Ga 



27r 

,-iHV'I + igs 



G 



R 



nKo{iR\/ri-iqs) 



G 



47tR 

The replacement 

cosh 29 



R 



4TrR 



Z,cos2e 



for ID 
for 2D 

for 3D. 

(72) 



followed by simple algebra makes the explicit expressions 
Ea. (j7ip for the generating function more compact 



P{C.V,z) 



dZ / d«;e-^^^^=(^^"-^) 



1 -1 
-2N'Ciqs(Z~w)(ri . ^_\2 



\Ga-GrY 



(73) 



At the limit z — > the position of the saddle point is 
independent from Z and w, which immediately yields 
the tunneling probability 



oo i 



and the number of the atoms 

P'AC.v,Q)= ^(Ga-Gr)' 



dZ I dwe-^'^'^^'^^-^Uz'-w') 



(74) 
(75) 



that have been kicked out from the condensate as the 
result of the tunneling. 

Calculation of the integrals Eq. ([71)1 for the probability 
results in 



^(C,^,o) 



{Ga — Gfl) ■ 



(A/'C?. 



2AA2C2(j2 ' 



(76) 



It is expedient to show in Fig llll the Fourier transforms 
over the frequency ^ of the coefhcients in front of the 
(Ga — Gr) 2 in Eq. ([7S|) . The electron tunnelling proba- 
bilitythrough the condensate linearly grows as function 
of time and saturates at the Heisenberg time, which is a 
natural consequence of the employed OD model of fluc- 
tuations (coordinate independent matrix Q). However, 
the typical saturation time t '^ Afqs depends on the en- 
ergy T] via the energy dependence of the position qs (77) of 
the saddle point position in the complex plane, which we 
will find later on. Note that, a rapidly oscillating part of 
the coordinate dependence, which is present due to the 
interference beats between Ga and Gr for the leads size 
small compared to l/y^, must disapear for large liads. 



w 



and the combination {Ga 
2GaGr. 



Gff) has to be replaced by 
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FIG. 11: Time dependence(not to scale), corresponding to the 
frequency dependence of the tunneling probability for 7 = 0. 



Let us now calculate the number of atoms kicked out 
of the condensate. At this stage we have to take into ac- 
count the fact that the frequency C, = uj — 'i'y has an imagi- 
nary part — i7 originating from the fact that the electron 
leaves the condensate once it has reached the outgoing 
lead. This implies that the electron spends just a finite 
period of time in the condensate and therefore kicks just 
a finite number of atoms out of the latter. Technically, 
this manifests itself in the divergency of the integral for 
P^(C, ?7, 0) emerging from Eq. (|75)) for the case 7 = 0. By 
direct integration over AZdw in Eg. ([75]) for the case of 
large leads one finds 



^;(C,^,o) = 



GaGr 



^2 {2MqsO' 



(77) 



Fourier transformation of Eas. (|76[l77l over ^ has the 
form 



P(i,7?,0) = GAGfl 



(ie (i) - (f - 2Uqs) e (t - 2Uqs)) 



2AA2g2 



P'At,v.Q)= - 






fQ [t] 



(78) 
(79) 



where 8 [x) is the step function. Note that the time 
dependence of the tunneling probability comes from the 
spectrum statistics; it is typical of the chosen Gaussian 
unitary ensemble (GUE). In Fig[T2]we present these two 
dependences for 7 7^ 0, The finite 7 results in the net 
tunneling probability 



7P(t,7/,0)df = GAGij 



1 _ e-^^T?' 
2Nqs 



and the total number of deplited atoms 



P^(i,r7,0)dMV 



Qs 



4 {M-iqs 



GaG RdV 



(80) 



(81) 



It corresponds to the situation where the destructive in- 
terference increases in the course of time and stops the 
tunneling process after a time of the order at the Heisen- 
berg time. At longer times, the electron does not present 
in the condensate and hence does not result in atomic 
losses. 



P, bMlMq, 




FIG. 12: Tunneling pronability P (solid line) and the number 
of particles 5N that have left the condensate in function of 
time line, (dashed line, not to scale) for 7 7^ 0. 



We now concentrate on the position of the saddle point. 
To this end one has to consider the saddle-point condi- 
tions suggested by the requirement 

da.Ws +/^ln[^-p'^-?.]U0, 



where the cut-off momentum is given by the condition 
,y ,d f^" d'^P = A/". One can consider the dimcnsion- 

(27r) J -^ 

less density n = J\f/V as a small parameter. Taking the 
derivative one finds 



2nqs 



dp' 



{2^t[^~p'^^q,) 



The left hand side of the equation contains the small 
parameter, and therefore the right hand side also should 
be small. The contribution to the integral on the right 
hand side comes from the domain of small momenta p' ^ 
n^l'^ while 77 ~ g^ '^ 1. We therefore can neglect p'^ in 
the denominator and arrive at 



2ngs 



dp' 



which immediately results in 



(82) 



i± vV^ 



(83) 



that is in the dependence typical of the random matrices 
and resulting in the Wigner semicircle law for the state 
density. Imaginary part of this quantity is the parameter 
igs entering Eas. (|72p . while the real part (which have 
been earlier ignored in the expressions for P and SAf) 
just gives a shift of the energy scale. Substitution of the 
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P(nM\ 




FIG. 13: Tunneling probability (bold lines) and the number 
of atoms kicked out of the condensate (thin lines) as functions 
of the energy rj (not to scale). Solid lines correspond to ID 
, dashed lines to 2D, and dash-dot lines to 3-D. Zero energy 
corresponds to the bottom of the band 77 = — -\/2. Here R = 10 
and 7 — 0.1. 



result Eq.dHSD to Eqs.d?!]) yields 



G 



A-R 



^2-1/2 



\]r\±^Jr^ 



: exp 



±ii? 



Ga-^h = ^K, I ±R [ -i#E! 



Ga;B. -^ 4^ exp 



±ii? 



-»j±V.?^-2 




for ID 
for 2D 

for 3D. 

(84) 



Now one has to substitute the obtained Green's func- 
tions to the expression Eqs. ljSQBip for the net tunnehng 
probabihty and the total number of atoms kicked out 
from the condensate. 



-fPit,7^,0)dt^GAGR 



1 



-2^1 qs 



2Nqs 



(85) 



and the total number of deplited atoms 



[ PUt,v,0)dtdV - f I GaGrAV 

J 4 [N-iqs) J 



(86) 



These results are shown in Fig[T3] One sees that the de- 
pendences -P(?y) are not too much sensitive to the dimen- 
sionality of the problem, and the main contributionto the 
conductance comes from the edge of the electron levels 
band, which is formed as the result of the perturbation 
of the mean-field band by the strong quantum fluctua- 
tions of the atomic density in the condensate. The de- 
pendences SAfiji) have more pronounced dependence on 
energy via dependence of the volume ocupied by the tun- 
neling electron on the spacial dimensionality. 
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